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Abstract 

We propose a tomographic scheme to reconstruct the quantum state of a Bose- 
Einstein condensate, exploiting the radiation field as a probe and considering the 
atomic internal degrees of freedom. The density matrix in the number state basis 
can be directly retrieved from the atom counting probabilities. 

PACS numbers: 03.65.Bz, 03.75.Fi, 32.80.-t 
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Since the Einstein's generalization of the black-body law derived by Bose, where the 
prediction that identical atoms considered as an ideal gas at sufficiently low temper- 
ature would become tight together in the lowest quantum state to form a condensate 
was made, there was no experimental evidence of this phenomenon, meanly due to the 
extremely low temperature at which it should occur. Finally, laser-based and magneto- 
optical trapping techniques were developped to cool and trap atoms at the microkelvin 
level, preventing their solidification and/or liquefaction. Then, in 1995 groups at JILA 
|1| and MIT obtained the experimental evidence for the condensation in dilute 
atomic gases: the Bose-Einstein condensation (BEC). The BEC is a macroscopic oc- 
cupation of the ground state of the gas and is one important paradigm of quantum 
statistical mechanics. The density distribution of the condensate should be represented 
by a macroscopic wave function with a well defined amplitude and phase, because of 
the transition from disordered to coherent matter waves. Very recently high-contrast 
interference of two independent condensates were observed showing that a conden- 
sate has a well defined phase. Although the evidence of the macroscopic phase seems 
well established by the experimental results of Ref. what is the BEC ground state 
is still an open question and one should devise a technique able to measure it. 

In the last years in quantum optics the state measurement of an electromagnetic 
field has become a very studied subject, both from theoretical and experimentral points 
of view, after the pioneering work by Vogel and Risken [Q. The fundamental idea is 
that in repeated experiments the outcomes of the homodyne measurement of the output 
light give the marginal distribution of the Wigner function of the field. By varying the 
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local oscillator phase one can reconstruct the whole Wigner function, then the state 
of the field. This is not, however, the only way to measure the quantum state of the 
radiation, because it was recently shown || ||, |7j that the reconstruction is possible 
by directly sampling the number of photons. Since the atoms in BEC are bosons as 
well, one can think to apply the same argument to reconstruct the condensate's state. 
Recently, another approach to reconstruct the BEC density matrix appeared in quant- 
phys files ||. The latter is however based on the knowledge of the relative phase of 
two non interacting condensates which are different for their hyperfine structures. It 
also needs of a beam splitter modeled for example by an output coupler for atoms P| . 
Instead, we shall consider a reconstruction scheme only based on atom counting, which 
exploits the interaction between the condensate atoms and a radiation field. 



To this end we follow Javanainen |T0[ and consider a Bose gas consisting of atoms 
moving in an isotropic harmonic oscillatory potential. The eigenstates of the center of 
mass motion |n) are labeled by the vector index n = (n x ,n y ,n z ) of positive integers. 
For the internal state of each atom we adopt the conventional two-state model: ground 
state \g), excited state |e) which we assume having a long enough life-time r a . These 
are separated by the frequency uj e ^ gi and the dipole moment matrix between the states 
is d. A radiation field, whose plane wave modes are enumerate by the index q, that 
incorporate both the wave vector q of the photon and the polarization vector e g , 
interacts with the ideal gas. The frequency of the mode q is fl q = |q| (we shall use 
here the convention % = c = 1). 

The Hamiltonian for the system consisting of the Bose gas and the photons thus is 
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n q 

~ E [e(g)(n'|e^ r |n)6t in ,6,, n a g + H.c.] . (1) 

The first two terms are the energies of the ground state and excited state atoms, the 
third term is the Hamiltonian for the free photon field and the final term governs the 
atom-field interactions. For instance, a process in which absorption of a photon in 
a state q converts a ground-state atom in the cm. state |n) to an excited atom in 
the cm. state |n'), is governed by the matrix element — £(g)(n'|e" 1 ' r |n). The coupling 
coefficient pertaining to the internal states is £(g) = ^Jfl q /2V e q ■ d, where V is the 
quantization volume. The explicit matrix elements in ([!]) is the overlap between the 
cm. state |n') and the cm. state |n) shifted by the momentum q. The matrix element 
d is taken such that the correct optical linewidth 7 = d 2 uj^ g /6ii results (for Cs atoms 
it is 7 « (2tt)2 MHz). 

In reality, for atoms moving inside a trap the total Hamiltonian should also contain 
terms representing the interaction between atoms, however, these can be neglected if 
we assume low enough density. We further assume that the condensate is probed by a 
weak laser beam whose mode index and frequency are denoted by k and Q. The latter 
will be resonant with the trapped atoms, i.e. Q = uj e ^ g ~ (2ir) 4.0 x 10 14 Hz for Cs 
atoms. As usual, the external field can be regarded as classical, and the corresponding 
photon operator can be treated as a c-number, i.e. we make the replacement —>■ a^. 
The spontaneous widths of the excited state are much larger than any relevant cm. 
frequency e n ; consequently we proceed from now on as if e n = 0. Then, the Hamiltonian 
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([]]) will become 

H = uj £ , 9 £ ^ An - Y: [ee-^in'le^bl^ + H.c] , (2) 

n n,n' 

where 9 = £,{k)\ak\ is the Rabi frequency that ensues when a classical field drives the 
internal transition in a single atom, and <fi = arga^. 

Usually the condensate is modeled by treating its annihiliation operator as a c- 
number, i.e. b g $ — > y/Nc with iV c the number of atoms in the condensate. Though this 
is true in the mean [[H], [11] , we continue to regard it as a purely quantum operator. On 
the other hand, as it results from the first term inside the square brackets of Eq. (Q), 
the driving field only couples states in pairs {\g) £g> |n) , |e) ® l^n}}, where |e) ® \tp n ) is 
obtained with the momentum translation k from the excited state |e) <g> |n). We now 
assume that (almost) all atoms are in the condensate, therefore all matrix elements 
referring to the state \g) ® |0) should pick up a large multiplier \/AQ. This allows 
us to retain in our theory only the condensate and its corresponding excited state 



e) ® |"0o) One then could resort to the boson operators 



b\ = b 9i0 (3) 
b\ = E<nV k ' r |0>&U (4) 

n' 

so that the mode 1 represents atoms in the internal ground state (and cm. ground 
state), while mode 2 represents atoms in the internal excited states (no matter of what 
cm. state). 

In view of the above assumptions and for Q = uj e ^ g) Eq. (||) can be simply written 
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as 



H = ttblh - 6 \blb ie - int+i4 + b 2 b\e tnt - i4 



(5) 



where the phase is governed by the applied field. 

If the atomic density operator at the initial observation time is p, then its evolution, 
in a frame rotating at the frequency of the applied field, will be given by 



p(r) = K(Q,4 > ypK(Q,4 > ) 



(6) 



with 



K(Q, <j)) = exp [ie {b\b 2 e i(t> + b^e* 



Q = 9r. 



(7) 



Soon after this evolution, the system is left to expand ballistically and the number 
of atoms are counted, with some selective detectors able to discriminate atoms in the 
excited state from those in the ground state. In order to simplify the presentation we 
assume that the measurement time is shorter than r a , otherwise the detection efficiency 
should be considered as in |L2|, 0] , while for sake of simplicity we assume unity detection 
efficiency. 

The evolution operator can be written in terms of generators of the SU(2) group 



J + = b\b 2 ; J_ = hb\ ; Jo = ~ (b\h - b\h 



(8) 



obtaining 



K(e, 



exp 



{'9 (j+e^ + J_e _<< 



(9) 



The probability of counting n\ atoms in the mode 1 and n 2 in the mode 2 at the 
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detector, for fixed d> and 0, is 



P(m,n 2 ,Q,(f)) = (n 1 ,n 2 |ir(9,0)V^(e,0)|n 1 ,n 2 ) 



E E K) (10) 

l\,m\ l2,m 2 



where 



^Kfa,™^) = (M 2 |^(©,0)Kn 2 )* x (m^m^Kie,^,^). (11) 

The rotation matrix elements in the number state basis are obtained by using the 
Baker-Campbell-Hausdorff formula for the Lie algebra of SU (2) p 



(m 1 ,m 2 \K( y Q,(f))\n u n2) = EE< 

j k 



\ (mi — fc)!m 2 !(rii — j)!n 2 ! 



X (Uanl " ; e)J+fc exp [- ln{cosh 0}(m - n 2 - 2j)l e^ fc "^ , 

where the prime sign on the sum means that k assumes the values n 2 — m 2 + j and 
m i — n i + j. 



Eqs. (11) and (12) show that 7Z is simply a linear operator connecting the density 
matrix elements to the number probability data which gives rise to an overdetermined 
system of linear equations provided one knows the probability as a function of the 



experimentally controlled parameters (p an d ©• Hence, it can be inverted to get 
the density matrix elements in the number basis by just collecting the experimental 
data for various values of and G. Since we do not assume any fixed total number 
of atoms, a truncation procedure should be employed. Let us label with N% and N 2 
the upper bounds for the number of atoms in the modes 1 and 2 respectively, then p 
should be intended as a complex vector with (Ni + l) 2 x (iV 2 + l) 2 elements and P 



as a real vector with (iVj + 1) x (N2 + 1) x x Nq elements, where and iVe are 
the total number of values of and G used. Provided to have > (JVi + 1) and 
> (-^2 + 1)) or viceversa, one can invert Eq. fllCf) to get the density matrix from the 
measured probabilities by using some numerical methods flSfl . Essentially, by defining 
the hermitian matrix G = TZ'TZ, we can compute 

p = G- 1 TZ^P, (13) 

provided to have G non singular. 

It is worth noting that the above introduced density matrix does not concern only 
the condensate atoms in the ground state, but also those in the excited states. Since, 
due to the low density assumption, we are neglecting the collisional effects between the 
two modes, it is reasonable to write it in a factorized form before the radiation field 
action, i.e. 

p = pM®pW, (14) 

where now p- l > represents the density operator for the condensate. Of course once one 
has calculated Eq. (|l"3"D, the matrix elements p^ mi immediately follow by tracing out 
the mode 2. 

The above procedure, however, although correct in principle, seems rather difficult 
to implement numerically due to the very high dimensions of the involved matrices. 
To simplify the scheme it is possible, with a good approximation, to assume that in 
the initial state all atoms are in their internal ground state, i.e. p^ = 1 0) 2 2 (0 1 , then 
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Eq. ( |TOD becomes 

p(m, n 2 , ©, 0) = J2 K;Zo,o( e > ^PhL ■ (15) 

In this case should be intended as a complex vector with (N% + l) 2 elements, then 
Eq. fll5|) can be inverted by only varying the phase <fi of the applied field, mantaining 
fixed, for example. The phase of the external field can be arbitrarily fixed with the 
initial measurement because the reconstruction procedure only needs to specify a given 
number of phases in the interval [0, 2n] starting from any phase value. 

In conlusion, we have presented an approach to the quantum state reconstruction 
of a BEC based on a tomographic procedure which seems easy to implement because 
one only needs to master the probing field and to measure the number probabilities 
of atomic detection. It is worth noting that the presented method does not need 
the knowledge of the exact number of the trapped number of atoms but only an upper 
bound for it. Of course, for a fixed and known number of atoms the procedure is simpler 
since one needs to only vary the phase of the field independently of the assumption on 
p&>. Only in this case, however, one can get rid of non unity efficiency of the detector, 
because by counting both species of atoms one can discard those data which do not 
give the initial total number ||. It is also worth noting that the presented procedure 
is well suited for condensates with a relatively small number of atoms, because of 
the number of data one should collect. This should not be a limitation because if 
the condensate were in an almost coherent state, as it is usually assumed, its wave 
function should not depend too much on the number of atoms in the condensate. 
Thus, BEC with a small number of atoms (say 50-100) should be obtained in order to 
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easily implement the present model, and in such a case the assumption we introduced 
might be better justified. The realization of such a small condensate could be the 
object of next generation experiments, once the existence of a given phase for the 
condensate will be better tested by other experimental groups. In our opinion the 
possibility of "measuring" the true density matrix of the condensate is accessible and 
worth considering. 

References 

[1] M. H. Anderson J. R. Ensher, M. R. Matthews, C. E. Wieneman, E. A. Cornell, 
Science 269, 198 (1995). 

[2] K. B. Davies, M. -O. Mewes, M. R. Andrews. N. J. van Druten, D. S. Durfee, D. 
M. Kurn, W. Ketterle, Phys. Rev. Lett. 75, 3969 (1995); C. C. Bradley, C. A. 
Sackett, J. J. Tollett and R. G. Hulet, Phys. Rev. Lett. 75, 1687 (1995). 

[3] M.R. Andrews, C.G. Townsend, H. J. Miesner, D. S. Durfee, D. M. Kurn, W. 
Ketterle, Science 275, 637 (1997), 

[4] K. Vogel and H. Risken, Phys. Rev. A 40, R2847 (1989). 

[5] S. Wallentowitz and W. Vogel, Phys. Rev. A 53, 4528 (1996). 

[6] K. Banaszek and K. Wodkiewicz, Phys. Rev. Lett. 76, 4344 (1996). 

[7] S. Mancini, P. Tombesi and V. I. Man'ko, Europhys. Lett. 37, 79 (1997). 

[8] E. L. Bolda, S. M. Tan and D. F. Walls, (|quant-ph/9703014j) 

10 



[9] M. -O. Mewes, M. R. Andrews, D. M. Kurn, D. S. Durfee, C. G. Townsend and 
W. Ketterle, Phys. Rev. Lett. 78, 582 (1997). 

[10] J. Javanainen, Phys. Rev. Lett. 72, 2375 (1994). 

[11] L. You, M. Lewenstein and J. Cooper, Phys. Rev. A 50, R3565 (1994). 
[12] U. Leonhardt and H. Paul, Phys. Rev. A 48, 4598 (1993). 

[13] Quantum Theory of Angular Momentum, L. C. Biedenharn and H. Van Dam, Eds. 
(Academic Press, New York, 1965). 

[14] D. R. Traux, Phys. Rev. D 31, 1988 (1985). 

[15] see e.g.: D. J. S. Robinson, A course in Linear Algebra with Applications (World 
Scientific, Singapore, 1991); W. T. Vetterling, S. A. Teukolsky, W. H. Press and B. 
P. Flannery, Numerical Recipes, (Cambridge University Press, Cambridge, 1985). 



11 



